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^ . ABSTRACT 

^v^j , We present 3D hydrodynamic simulations aimed at studying the dynamical and 

■ chemical evolution of the interstellar medium in dwarf spheroidal galaxies. This evo- 

\J3 ' lution is driven by the explosions of Type II and Type la supernovae, whose different 

OO . contribution is explicitly taken into account in our models. We compare our results 

CO ■ with detailed observations of the Draco galaxy. We assume star formation histories 

s ! | consisting of a number of instantaneous bursts separated by quiescent periods. Diverse 

, histories differ by the number of bursts, but all have the same total duration and give 

rise to the same amount of stars. Because of the large effectiveness of the radiative 
losses and the extended dark matter halo, no galactic wind develops, despite the to- 
tal energy released by the supernovae is much larger than the binding energy of the 
gas. This explains why the galaxy is able to form stars for a long period (> 3 Gyr), 
consistently with observations. In this picture, the end of the star formation and gas 
5^ ■ removal must result from external mechanisms, such as ram pressure and/or tidal in- 

^ | teraction with the Galaxy. The stellar [Fe/H] distributions found in our models match 

d • very well the observed ones. We find a mean value ([Fe/H]) = —1.65 with a spread 

of ~ 1.5 dex. The chemical properties of the stars derive by the different temporal 
evolution between Type la and Type II supernova rate, and by the different mixing 
of the metals produced by the two types of supernovae. We reproduce successfully the 
$H ' observed [0/Fe]-[Fe/H] diagram. However, our interpretation of this diagram differs 

from that generally adopted by previous chemical models. In fact, we find that the 
break observed in the diagram is not connected with the onset of a galactic wind or 
with a characteristic time scale for the sudden switchover of the Type la supernovae, 
as sometimes claimed. Instead, we find that the chemical properties of the stars de- 
rive, besides the different temporal evolution of the SNe II and SNe la rates, from the 
spatial inhomogeneous chemical enrichment due to the different dynamical behaviour 
between the remnants of the two types of supernovae. 

Key words: galaxies: dwarfs - galaxies: kinematics and dynamics - galaxies: abun- 
dances - galaxies: individual (Draco) - Local Group - hydrodynamics: numerical 
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1 INTRODUCTION 

Dwarf galaxies are t he most common class o f galaxies 
in the local Universe jMarzke fc da Costal Il997f) and they 
proba bly were much more numerous at past epochs lEllis] 
1997). Dwarf spheroidal/elliptical galaxies (dSph/dE) are 
found preferentially in high dens ity environments and can 
be stu died in nearby cluste rs dPopesso et alJ 12 0061) lik e 
Coma (ISecker fc Harrislll997l), Fornax jRakos et al.ll200lT . 
Virgo teinggeli fc Cameronlll99ll:lPhillipps et al.lll998l) and 
Perseus JConselice et al]l2003D . However, due to their prox- 



imity, dSphs of the Local Group (see Mateo 1998 for a re- 
view) offer an unique opportunity to study in detail their 
structural properties, formation and chemical evolution. 

Dwarf spheroidals are the least massive galaxies known, 
but yet, their velocity dispersions imply mass to light ra- 
tios as large as 100 Mq/Lq. This is usually explained as- 
suming that these systems are dark matter dominated. Ac- 
tually ,_jnth«_pastjew_y_ejxs_bj^ 

fe.g.lLokasll2002l: IWilkinson et al-lliool iKlevna et alJl2002l: 
Walk er et al .120051) and t h eoretical works(|K apmt^ Mjget^il] 
120041 iMaver et all l2005t iMashchenko et all l2005allbl) con- 
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firm the possibility that these galaxies are relatively mas- 
sive bounded system with virial masses in the rane e 10 s — 
5 x 10 9 Mp; b ut, fo r a dissonant view see iKrouoal l)l99?l) , 
iKroupa et alJ i2005l) . or the review by Gallagher & Wyse 
1994. 

Such galaxies are very metal poor and lack of neu- 
tral hydrogen and recent star formation. Thus they were 
initially believed to be very similar to Galactic globular 
clusters and to have a very simple star formation history 
(SFH). Recent studies have shown, instead, that these sys- 
tems are much more complex, with varied and extended 
SFHs. High resolution spectroscopy of severa l dSphs showed 
the p resence of a wide range in metallicity jHarbeck et alJ 
2001). For example, abundance analyses of stars belonging 
to Draco and Ursa Minor ha ve shown values of [Fe/Hl in 
the r a nge - 3 ^[Fe/H]^ -1.5 iStetsonlll98llShetrone et alJ 
Il998l . l200lT) with a mean value in the in terval —2.0 ^ 
([Fe/HD ^ —1.6, depending on the authors iAparicio et alJ 
1200 ll: IShetrone et~ai]|200ll : iBellazzini et al.ll2002l) ~ 

The above ranges are consistent, for some dSphs, with 
a single pe riod of star formation extended in time for 
a few Gvr dCarnev fe Seitzed Il986t iMateol Il998t iDolphinl 
| 2002l:lBabusiaux et al.ll2005h . As a further hint of long SFH 
Shct rone et al .1 teooifi found that their observed dSphs have 
[a/Fe] abundances that are ~ 0.2 dex lower than those of 
Galactic halo field stars in the same [Fe/H] range. This sug- 
gests that the stars in these systems were formed in gas 
pre-enriched by Type II supernovae (SNe II) as well as by 
Type la supernovae (SNe la), and star formation must thus 
continue over a relatively long timescale in order to allow 
a sufficient production of iron by SNe la (typically > 1 — 2 
Gyr). 

Given the small dynamical mass inferred for dSphs, the 
interstellar medium (ISM) binding energy is small when 
compared to the energy released by the SNe II explosions 
occurring during the star formation period; for instance, as 
shown in the next section, in a dSph as Draco, the baryonic 
matter has a binding energy of ~ 10 53 erg, while the ex- 
pected number of SN explosions in the past was 10 3 - 10 4 , 
realising an energy much larger than the binding energy. 
It is thus quite puzzling how the ISM can remain bound 
long enough to allow such a long star formation duration . 
Actually, cosmological simulations (e.g. lKawata et al"ll2005l : 
iRicotti fc Gnedinl2005l here after RG) find rather short du- 
rations (< 1 Gyr) of the star formation, and several au- 
thors believe that the ISM of these systems may be com- 
plete l y removed by violent SNe explosions fo ekel^^^iilk 



pletcly removed by violent blNe explosions (Uckel & but 
| l986|:|Mori et alll997tlMurakami fc Babull999HMori et aT 



2002, 200-tt lHensler et al.l2004ll together with the newly syn- 
thesized metals. 

Gas removal via galactic winds powered by SN 
explos ions is also invoked by lLanfranchi fc Matteuccil 
( 2004) (hereafter LM) who proposed a one zone chemical 
evolution model with very low star formation efficiencies and 
high wind efficiency, which is able to reproduce the metal- 
licity distribution and the general features of these galaxies. 
On the other hand, such features may be obtained also by 
closed box chemical evolution models by llkuta fc Arimotol 
( 2002) (here after I A) who prescribe a low star formation rate 
(SFR) and relatively long duration of the SFH, as suggested 
by observations. 

Motivated by the above arguments, in this paper we ex- 



plore the possibility that dSphs galaxies formed stars at a 
low SFR for a long period. To compare our results with ob- 
servations, we have tailored our models on the Draco galaxy: 
this galaxy is supposed to have experienced a star formation 
lasting for 3- 4 Gyr, and which essentially ceased 10 Gyr ago 
l)Mateoll998h . Obviously, our results may be confronted with 
other dSphs which are strongly dark matter dominated and 
have similar star formation histories as, e.g., Ursa Minor 
<Matedll99gl) . 

We run a number of three-dimensional (3D) hydrody- 
namical simulations to study the dynamical and chemical 
evolution of this system, following an assumed star forma- 
tion history. A special attention is paid to the influence of 
both SNe la and SNe II on the chemical enrichment of the 
new forming stars. 



2 THE MODEL 

In this section we describe as we build up our model. We 
start our simulations with the ISM in hydrostatic equilib- 
rium in the dark matter halo potential well. Any baryonic 
(stars or gas) contribution to the gravitational potential is 
neglected. We approximate the Draco observed stellar dis- 
tribution with a King profile in order to properly locate the 
forming stars and the SNe explosions on the numerical grid 
(see Section I2.2H . The dark halo properties are derived by 
the observed mass to light ratio, as described below. 



2.1 Galactic model 

We assume a quasi-isothermal dark halo, whose density is 
given by: 

Ph,0 



Ph(r) = 



1 + 



(1) 



where r^c is the core radius. The halo mass as a fuction of 
radius is then: 



M h (r) = 47rp h ,o?'h, c (y - arctant/), 



(2) 



where y — r/rh. c - The density profile is truncated at a tidal 
radius r^t in order to obtain a finite mass. The gravitational 
potential given by this mass profile is: 



$ h (r) = 47rGp hj0 r h ,, 



1 , .„ 2% arctany 
-log (1+y 2 + y - 

2 y 



(3) 



We assume the following relations among the dark halo 
parameters: 

1/3 



rh,, 



.9 x 10 



(i ' ■ V/| ' ^ 1 i, 

M 

-Mil \ , -2/3, 



h ' kpc, 



r M = 0.016 ( ^ ) /T^kpc 



(4) 
(5) 
(6) 



where h = 0.7 is the assumed value of the normalized 
Hubb le constant. The ab ove equations are obtained by 
iMac Low fc Ferraral lll999T) in the case of a quasi-isothermal 
model on the basis of similar relations found for dark matter 
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Table 1. Galaxy parameters 



Model 




Ph,0 


rh,c 


ni,t 


M 4 


P*,0 


r b 
cff 






(M/L v )g 


(M/L v )tot 




(10 6 M Q ) 


(1CT 24 g cm" 3 ) 


(pc) 


(pc) 


(10 6 M ) 


(10 -24 g cm- 3 ) 


(PC) 


(PC) 


(PC) 


(M /L v , ) 


(M Q /L V , ) 


Draco 


62 


4.3 


300 


1222 


0.56 


1.0 


210 


130 


650 


80 


211 


Draco S 


22 


(5.5 


160 


1005 


0.56 


1.0 


210 


130 


650 


14 


80 



All the above quantities are derived fro m the models, but: a the o bserved central stellar density iMateolll998t). b the observed 
effective radius of I he stellar component JPeterson fc Caldwelllll993ft . c the observed tidal radius Jlrwin fc Hatzidimit,rioiJll995h . d 



Value calculated inside the stellar volume. 



halos with Burkert p rofiles in th e range 4 x 10 9 < M h < 10 11 
M jBurkertj 1 19951) . F ollowing iMac Low fc Ferraral jl999l) 
and lSilich et alJ (l200ll) we scale down these relations to ha- 
los with lower masses. 

In the present paper we tailor our models on the Draco 
galaxy. The values of the parameters in the above formulae 
are constrained by observations of the mass to light ratio of 
this galaxy. We thus take into account the stellar component 
assuming a density profile given by the King profile: 



P*(r) = 



1 + 



3/2 ' 



(7) 



where r* iC is the stellar core radius. The stellar mass inside 
a given radius r is: 



(8) 



M*(r) = 47r/9*, rJ iC (\a{x - \j \ + x 2 ) 

where x = r/rv, c - As for the dark halo, the stellar profile 
must be truncated at a tidal radius r*,t in order to obtain 
a finite stellar mass. In the following we will refer to the 
volume inside r* it as the "stellar region", and to the volume 
inside rh,t as the "galactic region" . 

From observations of the central stellar l uminosity den - 
sity and the total V-band luminosity (e.g. iMated fl998l) . 
and adopting a stellar mass to light ratio in the V band 
M*/Ly = 2M /L , we estimate that p* : o = 1.0 x 
10~ 2 4 g cm" 3 a nd M* = M*(r*, t ) = 5.6 x 10 5 M 
(e.g. lMatedll998t) . We choose i %, c = 130 pc and r*^ = 
650 pc (see e.g. iLakd Il99(j IPrvor fc Kormendvl Il990t 
[ Irwin fc Hatzidimitrioul Il995t iPiatek et alJ 120011 ; IMated 
1998) which give the values of M* estimated above and 
also an effective radius r e g = 220 p c, in good agreement 
with the observed value r e s = 210 pc iPeterson fc Caldwell 



1993). Assuming that the total mass to light ratio in the 
volume up to r*, t is (M fe (r*, t ) +M*)/L V = 80 M /L v , 
jMatedll99sl) . we obtain the dark matter mass inside the 
stellar volume. In order to obtain the total dark halo mass 
Mh = Afh(fh,t) we solve iteratively equations 0] and 
[§] in order to find a dark halo with the right value of 
Mh(r*,t). Only one value of Mh exists, which gives the 
observed mass to light ratio. This value turns out to be 



halo masses 10 s — 


5 x 10 9 Mo (see e.g. [Klevna et al 


2002 


Kazantzidis et al. 


20, 0,4 1 Wilkinson et alJ|20Q4l: lLokas 


2002 


Walker et al. 2005; Mashchenko et al.ll2005al). 



We take into account the effect of the tidal interac- 
tion due to the Milky Way considering a somewhat less ex- 
ten ded dark halo than that indicated by equat ion El (sec 
e.g. iRead et alJl2005t iMaver et all 120051: iMashchenkoet alJ 



2005b). Thus we truncate the dark matter profile at rh,t = 
1.2 kpc (instead of 24.2 kpc) reducing the dark halo mass 
to a value of Mh = 6.2 x 10 7 M . This truncation does not 
affect the stellar component which is located in the inner 
region of the dark halo, a nd thus is much less affected by 
the tidal interaction (e.g. iGao et aT]l2004l) . The complete 
list of parameters concerning the galactic model is given in 
Table □ 

In this Table a model called Draco-S is also listed. This 
model has been considered to test the possibility that the 
dark halo could be lighter than that assumed above. Actu- 
ally, so me authors claim that dSphs could be n ot virialized 
(see e.g. lGallaeher fc Wvsell994l : lKrouDall997Fl . In this case 
the dark halo would be absent or with a mass similar to that 
of the stellar component. We did not consider such an ex- 
treme case but, in order to explore the effect of a lighter dark 
matter halo on our models, we run separated simulations in 
which M h = 2.2 x 10 7 M which gives a total M/Lv — 80 
M /L . As a consequence the dark halo mass distribution is 
more concentrated (cf. Table0 , while the stellar component 
is the same of the more massive models. 

We place the ISM in isothermal equilibrium within the 
potential well, with a temperature of Tism ~ T v i r - The ini- 
tial gas mass is Mism = 0.18Mh, which corre sponds to a 
baryonic fraction given bv lSpergel et alJ l)2006l) . The value 
of the central gas density pism.o is then fixed by the assump- 
tion that the ISM extends up to the dark matter component. 
The complete list of parameters concerning the ISM is given 
in Table H 

Finally, for computational reasons, we assume that the 
ISM has an initial metallicity Zjsm = 1O _4 Z . 



2.2 Supernovae explosions 

In order to understand the response of the ISM to the SFH 
and the SFR, we changed both of them in a number of mod- 
els all with the same initial conditions. We assumed that the 
star formation starts when the galaxy has accreted all its 
baryonic mass and that stars form in a sequence of A^burst 
instantaneous bursts, separated by quiescent periods. Sev- 
eral A^burst have been considered (see Table all giving rise 
to a final amount of stars equal to M* after f f = 3 Gyr, which 
is the end of the simulation. Obviously the lower is A^burst 
the higher is the intensity of each single burst (cf. Tabic 01. 
The prolonged star formation history adopted is consistent 
with re cent SFHs inferred for this galaxy by observations 
(see e.g. lCarnev fc SeitzeJlll8rilMatedll998l:lDolphinll2002h 
and other theoretical works (e.g TiltuTa^r^^motoll2002l) . 

After each starburst the stellar mass is increased by an 
amount AM* = M^/A^urst which is distributed in space fol- 
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Table 2. ISM parameters 



Model 


PISM,0 




M a 


M ISM,h 


bind 




(10- 24 g cm" 3 ) 


(K) 


(io 6 Mg) 


(io 6 M Q ) 


(10 51 erg) 


Draco 


0.4 


1.0 x 10 4 


3.3 


11.3 


83 


Draco S 


2.3 


2.0 x 10 3 


1.7 


3.8 


11 



a Amount of gas inside the stellar volume; b Amount of gas inside the dark 
halo volume; c Gas binding energy. 
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Figure 1. Logarithm of the total SN la rate T^g^ia for four dif- 
ferent star formation histories, all producing the same final mass 
of stars Mi, = 5.6 X 10 5 M after t f = 3 Gyr. Solid line: 50 
bursts; dotted line: 25 bursts; dot dashed line: 10 bursts; dashed 
line: single burst. 



lowing equation |7| We assume that the density profile of the 
stars does not change with time and that the newly form- 
ing stars are distributed according to the same nowadays 
shape; such an assumption is supported by se veral studies 
jMashchenko et aljl2005at iKawata et "al"ll2005l) . 

The SNe II explode at a constant rate 7?.snii for a period 
of t — 30 Myr (the lifetime of a 8 Mq star, the least mas- 
sive SN II progenitor) after each stellar burst. Thus the SN 
II rate for each burst is given by; 7?.snii = A?snii/t, where 
A^snii is the total number of SN II explosions occurring in a 
single burst. For a Salpeter initial mass function (IMF) it is 
iVsNii = 0.01 x (AM* /M Q ). Obviously AM*, M* and AWst 
satisfy the simple relation M* = A?burstAAi*. Our SFH also 
satisfies the relation U = NburstAt burst, where Afburst is the 
time interval between two successive bursts. Table [3] sum- 
marizes the models characterized by different choices of the 
above parameters. 

The SN la rate 7?.sNia for a si ngle population decreases 



in time after an initial rise fe.g. iGreggio fc Renzinil Il983t 
iMatteucci fc RecchilEoQll; iGreggid 1200511 . We adop t ed th 
time dependent rate given bv IMatteucci fc Recchil (12001 



according to the Single-Degenerate scenario (equation 2 of 
their paper with a fraction of binary systems A = 0.006). 
For every burst the time ti at which each SN la explodes is 
given by the condition: 
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Figure 2. Time evolution of the binding (solid line), thermal 
(thin solid line) and kinetic (dot dashed line) energies of the ISM. 
The upper panel refers to the reference model Draco-50, while the 
lower panel refers to the model Draco-25. 



where ti_i is the time at which the previous SN la (belonging 
to the same burst) exploded. 

In Fig. Qwe plot the time profile of the total SN la rate 
(given by the sum of all 7£sNia due to the bursts oc- 
curred within that time) for four different histories of star 
formation: a single burst, 10 bursts, 25 bursts and 50 bursts, 
all producing the same amount of total stellar mass M* after 
a time tf. For a single instantaneous burst the SNe la start 
occurring several Myr (> 30 Myr) after the onset of the star 
formation and continue for several Gyr owing to the long 
lifetimes of the stars responsible of this type of explosions. 
Note that when more than one burst is present TZg^ lBl has 
an oscillatory behaviour overimposed on the general increase 
with time. In contrast with 7£snii, 7?-sNia does not have an 
intermittent nature. This difference, together with the dif- 
ferent chemical composition of the ejecta of the two types 
of supernovae (see below), is responsible of some chemical 
properties of the stellar population (see Section [4.211 . 

Each SN explosion is stochastically placed into the 
galaxy. The probability of a SN explosion inside a radius 
r is proportional to the stellar mass within that radius: 
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P(r) — M+{r) / M+. The angular distribution is obtained 
choosing randomly the azimuthal angle $ and the cosine 
of the polar angle 8. 

Finally, we assume that each SN II ejects a mean mass 
of MsNii.cj = 10 Mq, and each SN la ejects MsNia.cj = 
1.4 M Q . Every SN II expels 1.0 M of oxygen and 0.07 
M[7) of iron (consistent w ith case A of IWooslev fe Weaver! 
Il99fil : iGibson et alJll99^ . and reference therein) while each 
SN la ejects 0.15 M n of oxygen and 0.74 M of iron 
iGibson et alj Il997l . and reference therein) . These yields 
slightly depend on metallicitv, un less metals are completely 
absent iWooslev fe Weaveil Il995h : we do not consider this 
dependence because it does n ot affect the [O/Fel feature s 
of the stars with [Fe/Hl > — 3 jOoswami fc Prantzosl boOO'). 
On the lGrevesse fc Sauvall jl998l) abundance scale, the iron 
and oxygen abundances for the SN II ejecta are Zo.snii = 
13.2 Zo,© and Zfg.snii = 5.7 Zf b ,q respectively, while for 
the SN la z7 ,SNia = 14.2 Z o ,0 and Z Fc , S Nia = 430 Z Fe ,Q. 
The explosion energy of each SN of both types is -Esn = 10 51 
erg. 



3 THE NUMERICAL METHOD 

We run the simulations with the 3D BOH (BOlogna Hy- 
drodynamics) hydro-code. The 3D code uses an Eulerian, 
second-order upwind scheme, in which consistent advection 
iNorman et al.lll980l) is implemented to reduce numerical 
diffusion. The 3D version of the code has been already used 
in different astrophysical co ntexts with satisfactory results 
jMarcolini et alJl2003Ll2004l) . 

In all the simulations we adopt 3D Cartesian coordi- 
nates. The total number of mesh points is 160 3 . The central 
volume is covered by 100 3 mesh points uniformly separated 
with Ax = Ay = Az = 13 pc; beyond this volume the linear 
mesh size increases geometrically in all directions with a size 
ratio of 1.14 between adjacent zones. In this way all the stel- 
lar region (where the SNe explode) is covered by the uniform 
grid and the grid edges are at large distances (~ 7 kpc) in 
order to avoid that possible spurious perturbations originat- 
ing at the boundaries may affect the solution in the central 
region. In all simulations outflow boundary conditions are 
enforced on all boundary planes. 

Besides solving the usual hydrodynamical equations, we 
also include two tracer variables representative of the SN II 
and SN la ejecta which are passively advected. Such tracers 
allow us to compute the chemical enrichment of the ISM due 
to the SN explosions. The ejecta (and the energy) of each 
SN explosion is initially distributed into a volume with a 
radius of two mesh points around the SN location. 

Radiative energy losses are taken into account consid- 
ering the metal dependent c ooling curve A as calculated by 
ISutherland fc Dopital (I1993T) . In any case the temperature 
is never allowed to decrease below 10 4 K. In the Draco- 
S models, due to the low virial temperature of the ISM, 
the gas is allowed to cool down to 10 2 K following, in 
this temperature range, the simple cooling curve given by 
iRosenberg fc Scheuerl lll973ft . 



4 THE REFERENCE MODEL: DRACO-50 
4.1 Hydrodynamical evolution 

It can be shown (see appendix A) that, given the SNe II 
rate for unit volume of our model, the supernova remnants 
(SNRs) overlap forming a network of tunnels filled by hot 
rarefied gas with large filling factors after nearly t — 12 Myr. 
Actually, our numerical model confirms that after this time a 
large fraction of the stellar volume is filled by the hot rarefied 
gas of the SNRs' interior while the dense SNR shells form 
dense cold filaments after colliding one with another. Despite 
the galactic gravity, the great majority of these filaments 
move outward pushed by the shock waves of the successive 
SN II explosions. In the merging cavities the gas has a mean 
density as low as 10 -28 g cm -3 and a temperature of few 
10 s K, and any further star formation is inhibited in these 
regions. Once the SNe II stop to explode (after 30 Myr since 
the beginning) the global cavity collapses and the ISM goes 
back into the potential well; this happens nearly 20 Myr after 
the last SN explosion. Note that the initial binding energy of 
the gas -Ebind ~ 8.3 x 10 52 erg is lower than the total energy 
1.12 x 10 53 erg released by the SNe II (see Tabled after a 
single burst. The simulation thus shows that the radiative 
losses are substantial and prevent the evacuation of the gas, 
as shown in Fig. [5] (upper panel) . In this figure the binding, 
thermal and kinetic energies of the gas are reported. It is 
apparent that these latter energies are much lower than the 
binding energy although the energy released by SNe II after 
only a single burst is 1.35-Ebmd- It is interesting to note 
that, despite its oscillatory behaviour, the mean amount of 
thermal energy of the gas remains constant, and does not 
increase in time as usually assumed in chemical evolutionary 
models. The consequences of this will be discussed further 
in Sections [U and [7| 

After 60 Myr the ISM distribution approximately recov- 
ers the initial condition, although turbulences and inhomo- 
geneities are now present (see Fig. [3] fourth column) . A sec- 
ond burst of star formation then occurs leading to a second 
sequence of SNII explosions. The gas undergoes a new cy- 
cle of merging bubbles which eventually collapse again. The 
influence of the SNIa explosions on the general hydrodynam- 
ical behaviour of the ISM is not very important because, as 
apparent in Fig. during a cycle of SN II explosions no 
more than 3-4 SNe la occur, only ~ 3% of the SN II num- 
ber. Despite their little importance from a dynamical point 
of view, the role of SNe la is very relevant for the chemical 
evolution of the stars (see next subsection). 

Figure |3] summarizes the evolution of the ISM and SN 
ejecta densities on the z — plane at several times. In the 
first and fourth rows the formation, expansion, and succes- 
sive interaction and re-collapse of the bubbles created by 
SNe II are illustrated. The central hole visible in the last 
panel of the first row is due to the explosion of a SN la oc- 
curred at t ~ 297 Myr. Analogously, the second and fifth 
rows show the evolution of the SN II ejecta density, while 
the third and sixth rows show the evolution of the SN la 
ejecta at the same times. From this picture it is apparent 
that the gas may leave the central region of the galaxy but 
it is never lost. At an early stage the ejecta is distributed 
quite clumpy, but becomes more and more homogeneous 
with time as more SNe II explode and the turbulence diffuse 
it through the ISM. Note that the SN la ejecta appears to 
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Figure 3. Logarithm of the density distribution (g cm -3 ) of the ISM (first and fourth rows), SN II ejecta (second and fifth rows) and 
SN la ejecta (third and sixth rows) in the z = plane at different times for the reference model Draco-50. The first, second, third and 
fourth columns represent snapshots of the gas after a time interval At = 15 Myr, 30 Myr, 45 Myr and 60 Myr from the occurrence of the 
latest instantaneous burst, respectively. Note that the time of the last column coincides with the occurrence of a new burst. Distances 
are eriven in knc. 
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Table 3. Supcrnovae parameters 



Model 


-^burst 


N a 

JV SNII 


(N S Nii/N burst ) b 


At£ t 

burst 

(Myr) 


N d 

IV SNIa 


Draco (S)-50 


50 


5.6 x 10 3 


112 


60 


254 


Draco (S)-25 


25 


5.6 x 10 3 


224 


120 


257 


Draco (S)-10 


10 


5.6 x 10 3 


560 


300 


266 



a Total number of SNe II; b Number of SNe II forming in each stellar burst; 
c Quiescent period between two consecutive stellar bursts; d Total number of 
SNe la after 3 Gyr. 
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Figure 4. Time evolution of the cold (T < 2 X 10 4 K) mass 
content inside the galactic (thick solid line) and stellar (thin solid 
line) regions. The dashed and dot-dashed lines refer to the warm 
(2 X 10 4 K < T < 10 6 K) and hot (T > 10 6 K) gas in the 
stellar region, respectively. The upper, middle and lower panels 
refer to the reference model Draco-50, Draco-25 and Draco-10, 
respectively. For the sake of simplicity only the evolution during 
the initial 300 Myr is reported. 



be distributed less homogeneously; the reason for this is the 
low SNe la rate. A deeper discussion on this point, and its 
consequences on the chemical evolution, will be presented in 
the next subsection. 

The periodic evolution of the ISM can be better seen in 
Fig. 2] (upper panel) where the behaviour of the gas content 
of the galaxy is plotted as a function of the time. In this 
figure three different phases of the gas are shown: cold gas 
with T < 2 x 10 4 K, warm gas with 2 x 10 4 K sC T ^ 
10 6 K and hot gas with T > 10 6 K. This distinction helps 
us to understand the thermal evolution of the gas and the 
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Figure 5. Time evolution of the mass content of the cold ISM 
(T < 2 X 10 4 K) inside the stellar (lower line) and galactic (upper 
line) regions for the reference model Draco-50 (upper panel) and 
model Draco-25 (lower panel). 



exchange among the cold, warm and hot phases. The warm 
and hot gas phases reside mainly in the stellar region, and 
their time evolution is indicated by the dashed and dot- 
dashed lines, respectively. 

For the sake of simplicity, in Fig. 2] we show only the 
evolution during five cycles in the time interval < t < 300 
Myr. The evolution over the entire period of 3 Gyr shows 
essentially the same periodicity, as can be seen in Fig. 
(upper panel) where the cold gas mass content of the galaxy 
is plotted up to 3 Gyr. 

Despite its oscillatory behaviour, nearly 30% of the 
initial gas is located beyond the galactic region. However, 
this gas is not lost by the galaxy but rather expands at 
larger distances remaining gravitationally bounded. Such 
a gas is obviously the most susceptible to be lost by ram 
pressure and tidal interaction with the Milky Way (see e.g . 
[ Marcolini et _alJl2003t ISofudll994l ; iMastropietro et aLllioOa ; 
iMaver et al.l2005ft . Indeed, since SNe feedback is not able to 
evacuate the ISM from the galaxy, in this scenario we must 
invoke an external mechanism, such as ram pressure strip- 
ping and/or tidal interaction, to remove the gas and thus 
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terminate the star formation. Ram pressure due to gaseous 
haloes of the Milky Way and M31 may also explain the ob- 
served correlation between stellar content and the gala cto- 
centric distance of dwarf galaxies llvan den Bergblll994l) . 

4.2 Chemical evolution 

The upper line in the top left panel of Fig. |S| represents 
the total amount of SN II ejecta expelled during 3 Gyr; the 
other two lines illustrate the evolution of the ejecta content 
inside the galactic (middle line) and the stellar (lower line) 
regions. The SN la ejecta has a similar behaviour (lower left 
panel). Note that along the entire evolution the fraction of 
the ejecta present inside the stellar region remains very low 
(~ 18% after 3 Gyr). This is the amount of metals which 
contributes to the metallicity of the forming stars. A large 
fraction of the ejecta is pushed at larger distances by the 
continuous action of the SN explosions. Indeed, nearly half 
of the ejecta still resides inside the galactic volume where it 
mixes with the ISM. The missing ~ 40% of the ejecta is still 
bounded and is mostly found in the neighbourhood of the 
galaxy, although a part of it is lost through breaches in the 
cold gas. 

Further insight on the chemical evolution is given by 
Fig. In the upper panel we plot the [Fe/H] distribution 
function (MDF) of the long-lived stars (with a mass ^ 0.9 
Mq), i.e. the mass fraction of these stars as a function of 
their [Fe/H]. The distribution is given at four different times 
in order to show its time evolution. 

We construct the MDF as follows. At every starburst 
the stellar mass in each ith computational cell of volume Vi 
is increased by an amount AM*,i = Ap*,iV!, where Ap*,i = 
P*,i/Mmrst represents the density of the newly formed stars 
(c.f. Section f2. 21 . Stars are not created in those cells where 
T > 2 x 10 4 K. The number of such cells in each starburst 
is however small and the general total stellar profile is still 
well described by equation [7| The newly formed stars are 
assumed to have the same chemical composition of the ISM 
in the ith cell at the moment of the burst: 

Z*,0 = (PISM^i°sm,/3 + PSNIlZsNII,/3 + psNIa^SNIa,/3)/Ptot , (9) 

where p to t = Pism + Psnii + psNia, and Z^ SM /3 is the initial 
ISM abundance of the generic element (5. 

At early times the MDF is rather broad reflecting an in- 
homogeneous distribution of the SN ejecta inside the galaxy 
because of the still low number of SNe explosions. As the 
number of the bursts increases, the ejecta becomes dis- 
tributed more and more homogeneously, and the MDF of the 
newly formed stars peaks around the mean value ([Fe/H]), 
which increases with time. For instance, the stars forming in 
the first burst are distributed rather shallowly in the range 
—4 <[Fe/H]< —1.9, while those formed in the last burst are 
clustered in the range —1.5 <[Fe/H]< —1.2. 

In Fig. |H] (upper panel) we show the final MDF together 
with the contribution due to SNe II only. Confronting these 
two types of contributions, we note that the effect of the 
enrichment due to SNe la produces a larger dispersion par- 
ticularly evident at the high [Fe/H] tail; this spread lowers 
the MDF maximum because the curve integrals (i.e. the to- 
tal mass of the stars) must be the same. 

The reason for the spread is given by the highly inho- 
mogeneous distribution of the SN la ejecta. Although the 



number of SNe la explosions after 3 Gyr is only 4.5% of 
that of SNe II, each of them expels 10 times more iron 
mass, and thus their contribution to the ISM metallicity 
is relevant. This contribution is particularly conspicuous in 
the high [Fe/H] tail of the MDF because, contrary to SNe 
II, SNe la may explode at any time, and in particular dur- 
ing the quiescent phases of SNe II. During these periods 
the ISM re-collapses and the SN la ejecta expelled in these 
phases remains more localized because of the larger density 
of the ambient medium. Actually, a rough calculation indi- 
cates that the SN la porosity (see appendix A) at this stage 
is Q ~ 0.04 (assuming a mean Z = 0.01 Zq). This means 
that the SNe la remnants are located quite apart one from 
another. Stars forming in these regions enriched by SNe la 
have an enhanced [Fe/H] ratio compared to the stars form- 
ing in the rest of the stellar volume, and are responsible of 
the high [Fe/H] tail of the MDF. 

At the end of the simulation we obtain a mean value 
of ([Fe/H]}=-1.65 with a spread of ~ 1.5 dex, in rea- 
sonable agreement with observations |Le hnert et al-|[l992|; 
IShetrone etafl 1200 it lAparicio et"al"ll200lt I Bclla zzini et alJ 
2002), while the maximum value of the distribution occurs 
at [Fe/H] ~- 1.5. 

In the lower panel of Fig. HQ t he final MDF is co mpared 
with the MDF obtained b y |B^UBzzini_et_aQ l|2002t) , u sing 
both the abundance scale bv lCarretta fe Grattonl dl997f) and 
IZinnfc WeslJl984l) . The calculated MDF is in good agree- 
ment with the observed one at the high [Fe/H] end. The 
tail of our MDF in the low [Fe/H] range is instead higher 
and less steep than the observed MDF. This is likely due to 
our simple SFH (bursts of equal strength) assumed. A more 
realistic SFH, with weaker early bursts, could depress our 
final MDF at low [Fe/ H], Also, the rela t ive in sensitivity of 
the method adopted bv lBellazzini et all i2002T) to the most 
metal poor stars, possibly reduces the discrepancy indicated 
in Fig. H 

ThE effect of the MDF tail at high [Fe/H] discussed 
above is particularly evident in Fig. [5] where the diagram 
[0/Fe]-[Fe/H] is shown at four different times. The open cir- 
cles form a statistically representative sample of the stellar 
distributions in [O/Fe] and [Fe/H]. At each time the sample 
is composed by a fraction 10~ 4 of the total number of stars. 
The solid line represents the mean value of the [O/Fe] distri- 
bution for any fixed [Fe/H]. The dotted lines in the fourth 
panel indicate the 1 a spread of this distribution. At any 
time a plateau at low values of [Fe/H] is delineated by the 
solid line, followed by a sharp decrease. With time the break 
in the distribution moves toward higher values of [Fe/H]. 
The plateau is due to the initially dominant contribution 
of SNe II to the Fe abundance: being the oxygen formed 
mainly by SNe II, the ratio [O/Fe] must remain nearly con- 
stant as the iron increases. The small negative gradient of 
the plateau is due to the slowly growing contribution in the 
Fe enrichment by SNe la. The sharply decreasing branch 
at higher [Fe/H] is due to stars formed in the regions of 
ISM recently polluted (mostly by iron) by SNe la. For this 
reason the stars populating the solid line in the decreasing 
branch are distributed diagonally nearly following a -1 slope. 
A glance at Fig. |7| and Fig. [5] shows that the stars on the 
decreasing branch populate the MDF high [Fe/H] tail, while 
the majority of the stars occupies the high [Fe/H] edge of 
the plateau in the [0/Fe]-[Fe/H] diagram. 
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Figure 6. Time evolution of the mass of the SN II (upper panels) and SN la (lower panels) ejecta inside the stellar (lower lines) and 
galactic (middle lines) regions. The upper stair shaped lines in all panels represent the total amount of SN ejecta expelled during 3 Gyr. 
The first and the second columns refer to the models Draco-50 and Draco-25, respectively. 



In the fourth panel of Fig. data from IShetrone et alJ 
teOOll) are reported. The four stars seem more uniformly 
distributed in [Fe/H] than our sample. This is because these 
four stars have been selected in order to span the largest 
range in [Fe/H] (Shetrone, private communication); thus 
their distribution has no statistical meaning and is not ex- 
pected to follow our sample distribution. We note that the 
point with the lowest [O/Fe] is displaced by ~ 0.5 dex 
to the left relative to the solid line denoting the descend- 
ing branch. This discrepancy does not represent a severe 
drawback for our model; there are in fact several explana- 
tions for it. On the theoretical side, had we started with a 
more extended gaseous halo, we would have found a lower 
mean value of [Fe/H], moving our sample to the left in 
the [0/Fe]-[Fe/H] diagram. On the observative side, the de- 
rived mean [Fe/H] value depends on the metallicity scale 
adopted (cf. the lower panel of Fig. 13 : for Draco it range s 
in the interval -2 < ( [Fe/H]) < -1.6 llAparicio et aljl200ll) . 
IShetrone et alJ i200ll) (whose data reported in our figure 
are, to our knowledge, the only published for Draco) find 



([Fe/H]) = —2.0: it is reasonable to hypothesize that, with 
a different choice of the metallicity scale, the observed stars 
reported in Fig. [5] would shift to the right by an amount 
very nicely fitted by our model. 

It is interesting to analyse the radial distribution of 
the ratio [Fe/H] within our model. Age/me tallicity gradi- 
ents h ave been c l aimed by some authors jHarbeck et alJ 
l200lt iLee et alJ 12001 iBellazzini et alJ 120051) to be 
present in several dSphs while others do no t find 
them dHurley-Keller et alJ Il999t Aparicio et alJ 1200 ll : 
ICarrera et alJl2002l : iBabusiaux et alJl2005l) The radial de- 
pendence of the ratios [Fe/H] and [O/H] at three different 
times are shown in Fig. 1101 where the circles form a statis- 
tically representative sample of the stars and the solid lines 
represent the radial profiles of the mean values of [Fe/H] and 
[O/H]. A big metallicity spread of ~ 1.5 dex is pr esent at 
any r adius (consistent with some observations, e.g. iTolstovl 
2004), and the profiles have been obtained splitting the 
galaxy in 50 spherical shells and averaging inside each shell. 

A slight metallicity gradient is actually present and re- 
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Figure 7. [Fe/H] distribution function of the long-lived stars formed up to four different times for the reference model Draco-50 (upper 
panel) and model Draco-25 (lower panel). Dotted line: 240 Myr; dashed line: 600 Myr; dot dashed line: 1.5 Gyr; solid line: 3 Gyr. 



duces with time. After 3 Gyr the difference between the cen- 
tral and the o uter region is 0.2 dex in [Fe/H], consistent with 
observations llHarbeck et all200ll) , while the [O /H] gradient 
is virtually absent. The origin and the evolution of the gra- 
dients are given by the following reasons. The ejecta are 
subject to two opposite effects: on one hand they are re- 
leased from SNe distributed following a King profile, thus 
more concentrated toward the centre, and on the other hand 
the hydrodynamics tend to distribute them more uniformly. 
With time the mean gas metallicity increases while becom- 
ing more and more uniformly distributed in space, and the 



metallicity gradient, essentially given by the freshly expelled 
metals, becomes less and less important. 



We point out that the [O/H] gradient is shallower than 
that of [Fe/H]. In fact, a substantial amount of iron present 
in the ISM is produced by SNe la whose ejecta, as discussed 
above, tends to remain more localized toward the centre 
(in particular during the quiescent phases) than the oxygen, 
which is expelled mainly by SNe II. 
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Figure 9. Abundance ratio [O/Fc] plotted against [Fe/H] of jV s sampled stars at four different times for the reference model Draco-50 
(left column) and model Draco-25 (right column). At each time N a represents a fraction 10 — 4 of the total number of stars at that time. 
The solid line represents the mean value of the [O /Fe] distribution for any fixed [Fe /H] . T he dotted lines in the f ourth panels represent the 
1 <T spread of this distribution. In the fourth panels we also show the values obtained by Shct rone et alJ l200llK In the lowest panels the 
same diagrams are obtained mimicking a one- zone model, i.e. assuming that the stars forming at each time have all the same metallicity, 
given by the spatial mean metallicity of the ISM at that time. 



5 OTHER MODELS 
5.1 Draco-25 

In order to understand the influence of the SFH on the chem- 
ical and dynamical evolution of the ISM, we computed a 
model similar to the reference one but with a time interval 
Atburst = 120 Myr between two successive bursts. As a con- 



sequence 7£snii = 7.47 x 10 6 yr 1 , doubled respect to the 
rate of the reference model. 

Due to this higher rate, the size of the cavities forming 
in the ISM inside the stellar region are larger (see Fig. 1111 . 
The gas at the periphery of the galaxy results to be more 
extended and less dense than in the reference model (com- 
pare Fig. 1111 and Fig. However, although the total SN 
energy released during one burst is 2.7 times the initial ISM 
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Figure 10. Radial distribution of [Fe/H] (left panels) and [O/H] (right panels) at three different times for the reference model Draco-50. 
The solid lines represent the radial profile of the mean values. 7V S has the same meaning as in Fig. |5] 



binding energy, radiative losses are still effective enough to 
impede the gas to be lost by the galaxy. Given the larger 
value of Atburst, the ISM has more time to pile up in the 
central region during the re-collapse before the next burst, 
thus reaching higher densities than in the Draco-50 model. 
In this region, which has actually the dimension of the stel- 
lar region, conditions are recovered for a new burst of star 
formation. 

The evolution of the ISM mass content inside the galac- 
tic and stellar regions is shown in Fig. |S] (lower panel) . Con- 
trary to the reference model, here the oscillations are much 
more regular because of the larger value of Atb urs t- 

The fraction of the ejecta present inside the stellar re- 
gion at the end of the simulation (see Fig |SJ is slightly 
smaller than in the reference model. Indeed, being the sin- 
gle bursts more powerful, the hot metal rich gas has more 
chances to leave the galaxy. As a consequence, the MDF 
peaks at a lower value of [Fe/H]~ —1.6. 

The chemical enrichment of the galaxy is synthesized 
in the lower panels of Fig. and Fig. |H1 and the right col- 
umn of Fig [5] We obtain a final value ([Fe/H]) — 1.7 with a 
spread of ~ 1.5 dex. The descending branch in the [O/Fe]- 
[Fe/H] diagram is now closer to the star with the lowest 
value of [O/Fe]. In general, however, these results show that 
the chemical evolution is slightly dependent on the adopted 
SFH. 

In Fig. 1121 we plot the final radial distribution of the 



ratio [Fe/H] and [O/H], as well as the radial profile of their 
mean values. As for the reference model, the chemical gra- 
dient results to be very shallow at the beginning of the sim- 
ulation, and tends to disappear with time. 



5.2 Draco-10 

In this model only ten bursts occur during 3 Gyr. They are 
separated in time by 300 Myr, and each of them is five times 
more powerful than those of the reference model. 

Each single burst tends to expand the ISM more than 
in the previous models, but radiative losses are still able to 
impede the gas to leave the galaxy. Indeed, as apparent in 
Fig. [I] (lower panel) , the mass content of the cold gas inside 
the stellar region may become as low as 10 4 M© just after 
the SNe II activity, but then grows again to the initial value 
during the SNe II quiescent phase. The gas continues to pile 
up toward the centre until, at t ~ 150 Myr, rebounds and its 
content decreases slightly (see the depression of the temporal 
mass profile at t = 250 Myr). 

After 600 Myr the chemical properties of the forming 
stars differ only slightly from those of the previous mod- 
els at the same time, indicating that the chemical evolu- 
tion is rather insensitive to the chosen SFH. For this reason, 
given the huge amount of CPU time required to complete 
the model, we stopped our simulation after t = 600 Myr. 
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Figure 8. Upper panel: [Fe/H] distribution function of the stars 
at the end of the simulation for the reference model Draco-50 
(solid line). The dashed line illustrates the contribution due to 
SNe II only. Lower panel: normalized MDF at t=3 Gyr for the 
models Draco-50 (thick solid line) and Draco-25 (thin solid line) 
comp ared to the observative MDFs of Draco (Bcllazzini ct al. 
2002) in the Zinn-West (dashed line) and Carretta-Gratton (dot- 
ted line) metallicity scales. 



5.3 Draco-S 

As discussed in Section |5J these models have been run 
to test the possibility that the dark halo could be lighter 
than the one we assumed in the reference model. In 
fact, relaxing the hypotheses of isotropic velocity dis- 
persions and virial equilibrium, one obtains mass-to-light 
ratios much lower than commonly ado pted (see e.g. 
iGomez-Flechoso fc Martfnez-Delgadol 12003. for the case of 
UMi). The non baryonic mass content (and consequently 
the baryonic one) in the model Draco-S is assumed to be 
only one third of the st andard mod el (c.f. Table 1 and 
2): M h = 2.2 x 10 7 M Q (lMatedll99l) which gives a total 
M/L v = 80 M /L Q . 

We explored the same SFH already discussed for the 
previous models with iVburst — 10, 25 and 50. In Fig. I13l the 
evolution of the mass content inside the dark matter halo is 
shown for all the three models. The galaxy gets rid of its ISM 
in all models on a timescale which increases with the number 
of bursts. This is due to the fact that when only few bursts 



occur, each of them is rather powerful (c.f. section 2.2), and 
can eject the ISM sooner than many weaker bursts. In all 
cases the gas is lost before t ~ 250 Myr, a time too short 
to be compatible with th e star format i on duration in ferred 
by the observations (e.g. lMateclll998t lDolphinll2002l) . Our 
results (in the limit of their assumptions, e.g. the ISM mass 
content and distribution are related to the DM ones and self 
gravity is neglected) indicate that the total M/L must be 
quite high and the dark halo must extend far beyond the 
stellar component. 



6 COMPARISON TO OTHER MODELS 

The SFH and the chemical enrichment of Draco and other 
dSphs have been discussed in terms of purely one zone chem- 
ical or cosmological models by several authors. 

• In their one zone model LM find that the plateau in the 
[0/Fe]-[Fe/H] plane is representative of SNe II enrichment, 
although their plateau level is higher than that found in our 
models because LM are able to distinguish the contribution 
to the metal production by SNe II of different masses, while 
we assume that every SN II has the same (average) chemical 
composition. The sudden [O/Fe] decrease at larger values of 
[Fe/H] is interpreted by LM as a consequence of a galactic 
wind which decreases the amount of available gas; the for- 
mation of new stars is thus slowed down and so is the pro- 
duction and the injection of a elements into the ISM by SNe 
II. SNe la, instead, h ave progenitors with much longer evo- 
lutionary timescales iMatteucci fc Greggiolll986l) . and con- 
tinue to inject iron lowering the [O/Fe] ratio. On the con- 
trary, in our models we do not have the formation of any 
wind and the change of the slope in the [O/Fe] vs. [Fe/H] 
diagram is given, as discussed in Section 14.21 by the inho- 
mogeneous distribution of the SNIa ejecta. In both models 
the descending branch is less populated than the plateau. 
In the LM models this happens because the most metallic 
stars form during the wind phase, when the star formation 
rate is greatly reduced; more metallic stars can be obtained 
increasing the star formation efficiency, although this leads 
to an overestimate of the [O/Fe] ratio. 

It is interesting to discuss further the differences between 
our models and those of LM. In these latter models the star 
formation rate depends on the ISM density and varies with 
time, while the SN II efficiency in heating the ISM has a 
fixed value assumed "a priori" (esnii = 0.03 in their no- 
tation), i.e. only a small fraction esnii of the energy of the 
SNe II explosions is transferred to the surrounding gas, the 
rest being radiated away. This energy fraction is constantly 
accumulated in the ISM, and the occurrence of a galactic 
wind at a certain time is unavoidable. The choice of the 
value of esnii is crucial for the behaviour of the models. 
An increase/decrease in its value anticipates/retards the oc- 
currence of the wind (for a fixed potential well). Thus, an 
appropriate choice of the ratio between thermalized SN II 
energy and potential well depth is required to obtain a rea- 
sonable duration of the star formation. Our models, in a 
sense, work opposite to those of LM. Intensity and duration 
of the star formation are fixed "a priori" , while the fraction 
of SNe II energy transmitted into the ISM is regulated by 
hydrodynamics and radiative losses. In our models winds ei- 
ther occur very soon or do not occur at all. In this scenario, 
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Figure 11. Logarithm of the density distribution (g cm -3 ) of the ISM (first and fourth rows), SN II ejecta (second and fifth rows) and 
SN la ejecta (third and sixth rows) in the z = plane at different times for the model Draco-25. The first, second, third and fourth 
columns represent snapshots of the gas after a time interval Af = 15 Myr, 30 Myr, 50 Myr and 80 Myr from the occurrence of the latest 
instantaneous burst. Distances are given in kpc. 
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a prolonged star formation period as that inferred for Draco 
(and other dSphs) suggests that the termination of star for- 
mation is due to an external cause such as gas stripping by 
the Galaxy (e.g. iMaver et allfeoOSl) . 

• One zone chemical models for Draco (and other dSphs) 
have been presented also by IA who do not take into account 
the possibility of galactic winds and assume that the gas re- 
moval should eventually result from an external mechanism, 
just as in our scenario. In absence of a wind, the IA interpre- 
tation of the sudden decline of [a/Fe] in the [a/Fe]-[Fe/H] 
diagram is based on the different temporal behaviours of the 
SN II and SN la rates. Given the much longer evolutionary 
timescales of these latter, an (assumed) sudden switchover 
of the iron source from SNe II to SNe la creates the break 
in the [a/Fe] ratio at a certain [Fe/H] value, corresponding 
to the time ii a at which a significant number of SNe la is 
assum ed to explode . This time was evaluated to be £i a ~ 1.5 
Gyr bv lYoshii et al.l lll99ft) who studied the stars in the solar 
neighbourhood. Before ti & the effect of SNe la on the metal- 
licity is not considered; for this reason the plateau obtained 
by IA is totally flat (see their Fig. 1). After 1.5 Gyr the Fe 
ejected by SNe la produces a decline in [a /Fe]. However, 
this decline is less pronounced than in our models and in 
those of LM. Indeed, such a decrease must not be confused 
with the decreasing branches obtained by the models of LM 
and ours. Rather, it is equivalent to the slow decrease of the 
plateau shown in these latter models where the action of the 
SNIa starts (and grows) just after 30 Myr. As a check, we 
calculated the average ratios [O/Fe] and [Fe/H] inside the 
stellar volume mimicking an one zone model. The results are 
shown in the lowest panels of Fig. |2l no break is present and 
the plateau has a slope similar to the mean slope of IA. In 
order to take into account the larger decrement of [a/Fe] ob- 
served IA must invoke a stellar mass function steeper than 
the Salpeter one, an assumption not necessary in our mod- 
els. 

• Cosmological simulations of the dynamical and chemi- 
cal evolution of dSphs have been recently performed by RG. 
These authors point out that objects with masses below 
~ 10 s M are able to form stars before the reionization 
epoch (12.5-13 Gyr ago) losing their gas because of the pho- 
toheating by massive stars inside the galaxy (and not by 
the external ionizing background). More massive objects, 
instead, retain a fraction of their gas even after the reion- 
ization. Thus, the duration of the star formation in small 
objects is rather short (< 1 Gyr), while in more massive 
galaxies it can be protracted for a longer time. RG iden- 
tify these latter galaxies with the observed dwarf irregulars, 
while the observed dSphs properties are reproduced by the 
smaller simulated objects. These latter objects show a metal- 
licity spread that can be understood in terms of hierarchical 
accretion of subhalos containing stars with different metal- 
licities. RG thus argue that the observed metallicity spreads 
are not necessarily an indicatio n of star formation extende d 
over several Gyr, as claimed bv lGrebel fc Gallagher! l)2004l) . 

Comparing the above results with ours, we point out that 
the longer star formation assumed in our models (with Mh = 
7 x 10 8 Mq) is marginally consistent with the finding of RG. 
We also stress that, although the metallicity spreads found 
in our models derive from a completely different mechanism, 
yet we confirm that large spreads of [Fe/H] can rise quite 
early in our simulations (cf. Fig. 7). 
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Figure 13. Time evolution of the mass of the cold ISM (T < 
2 X 10 4 K) inside the galactic region for the Draco-S models. The 
mass is normalized to the initial mass inside the galactic region. 
Solid line: Draco-50; dashed line: Draco-25; dotted line: Draco-10. 
The stars indicate the time at which the stellar bursts occur. 
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Figure 14. Central distribution on the z = plane of the SN 
II (upper panel) and SN la (lower panel) ejecta densities for the 
Draco-25 model at t = 1520 Myr. These panels are a close up view 
of the last two panels of the fourth column in Fig llll Distances 
are given in kpc. 
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Figure 12. Radial distribution of [Fe/H] (left panel) and [O/H] (right panel) at the end of the simulation for the model Draco-25. The 
solid lines represent the radial profile of the mean values. N s has the same meaning as in Fig. |3 



On the other hand it is known dShetrone et aljfeoO lT) that 
dSphs have [cc/Fe] ratios ~ 0.2 dex lower than those of the 
Galactic halo fields s tars in the same [Fe/ Hl range, whose 
[a/Fe] ratios are 0.3 dShetrone et al.M2001f) . This usually is 
interpreted as a signature of a longer star formation occurred 



SNe la jMatteucci 


2003). Alternative suggestions based on 


different IMF (e.g. 


Shetrone et al.||200l|) or on less rotation 



of massive stars ( compared to sol ar neighbourhood) have 
also been invoked dTsuiimotoll2006l) . but we do not consider 
here these speculative hypotheses. 

A glance at Fig. |U] illustrates that in our simulations a 
decrement of ~ 0.2 dex in the plateau is achieved after 
t ~ 2.0 Gyr. Thus we conclude that, at least in our models, 
only a prolonged star formation history (> 2.0 Gyr) can ac- 
count for the chemical differences between the Galactic halo 
and the dSphs. Such a long SFH for dSphs is also found in 
the chemical models bv lFenner et all fcOOfj) who are unable 
to reproduce the Ba/Y ratio unless stars formed over an in- 
terval long enough for the low-mass stars to pollute the ISM 
with s-elements. 

Fig. [5] also shows that the formation at early times of 
stars with low [O/Fe] and relatively large [Fe/H] is not pre- 
cluded, although must be regarded as a rather rare event. 
Thus our models do not exclude short star formation du- 
rations, although the formation of stars with low [a /Fe] is 
less probable. 

• Cosmo logical simulations o f dSphs have been worked 
out also bv lKawata et al] <l2005l) who, in line with RG, find 
a rather shor t (< 1 Gyr) duration of the star formation. 
iKawata et alJ <l2005t) focus on Sculptor galaxy in which Tol- 
stoy et al. (2004) found two distinct populations, with the 
lower metallicity stars more spatially exte nded than the 
higher metallicity stars. This is interpreted bv lKawata et alJ 
(2005) as a consequence of the metallicity gradient which is 
present in their simulated system. This gradient is similar to 
those found in our simulations at early times; in our mod- 
els, however, the gradients decrease with time. We thus find 
that, if two distinct stellar populations arise from a metal- 
licity gradient (different possible mechanisms leading to dif- 
ferent populations are given by Tolstoy et al. 2004), then 
their presence is indicative of a short star formation history. 

We conclude ou r comparison with the model by 
IKawata et alJ <l2005l) briefly highlighting two differences: i) 
ou r models do n o t get rid of the gas, while the in model 
bv IKawata et all (120051) the SN feedback leads to a blow- 
out phase; u) no break is present in the [0/Fe]-[Fe/H] dia- 



gram obtained bv lKawata et all (|2005). The first difference 
is due to the initial amount of gas assumed by these authors, 
which is an order of magnitude lower than ours; as a conse- 
quence, the gas binding energy is much lower and the radia- 
tive losses (which are proportional to the square of the gas 
density) are greatly reduced, allowing the SNe to vent the 
gas away. The secon d difference comes, as fully realized by 
IKawata et alJ $2005), from having t hese authors assumed for 
the SNe la the model proposed bv lKobavashi et al.l <l2000l) 
who suggested that SNe la are inhibited in the stars with 
[Fe/H]< -1. 



7 DISCUSSION AND CONCLUSION 
7.1 Dynamical evolution 

We run several 3D numerical models in order to understand 
the chemical and dynamical evolution of the ISM in dSphs 
as a function of several parameters like the SFH and the 
amount of dark matter. Given the large amount of cpu time 
required to complete a single model, we could not explore 
a vast region of the parameter space; instead, we tailored 
our galaxy model on Draco and considered SFHs similar to 
those supported by the observations for this galaxy. All the 
adopted SFHs last 3 Gyr and give rise to the same final 
amount M* — 5.6 x 10 5 Mq of stellar mass. The stars are 
supposed to form in a sequence of instantaneous identical 
bursts separated by quiescent periods; different SFHs differ 
by the assumed number of bursts, but the total number of 
SNII (but not SNIa) is always the same. 

Although the total energy released by the SNe is orders 
of magnitude larger than the binding energy of the ISM, the 
galaxy retains almost the totality of its initial gas, unless 
rather light dark haloes are assumed. This is due to the 
huge efficiency of the radiative cooling, despite the low mean 
metallicity [Z ~ 10~ 2 Zq) of the gas. For most parameter 
combinations tested in this study, such effective radiative 
losses prevent the gas to slowly accumulate the SN energy. 
With our assumptions, the galaxy never gets rid of its gas, 
unless the potential well is rather shallow (see Section f5.3H : 
in this latter case a galactic wind starts quite soon and the 
galaxy looses all its gas in less than 200 Myr. 

In the light of the above discussion, we conclude that 
the dark halo of Draco must be massive and extended in 
order to retain the gas for a period of several Gyr, the du- 
ration of its star formation. This in turn implies the need of 
an external mechanism to remove the gas and end the star 



Star formation feedback and metal enrichment in dSphs 17 



formatio n, as gas stripping and/or tidal interaction by the 
Galaxy (iMaver et alJl2005l) . 

7.2 Radiative losses 

The above arguments highlight the crucial role played by 
the radiative cooling in our models: would these losses be 
less substantial, galactic winds could develop more easily 
and the evolutionary scenario would drastically change. A 
correct evaluation of the radiative losses in the numerical 
simulations is thus of the utmost importance. Unfortunately, 
it is known that in numerical hydrodynamics several factors 
concur to degrade an accurate estimate of these losses. More 
confusing, some factors lead to an overestimate, while others 
to an underestimate. 

Overestimates occur particularly at the contact discon- 
tinuities separating hot rarefied and cold dense gas phases 
because the intrinsic diffusion of the numerical scheme 
spreads these discontinuities over several mesh points cre- 
ating a gas phase with intermediate densities and tempera- 
tures characterized by a large emissivity. Although several 
physical processes (such as turbulences and heat conduc- 
tion) really smear out these discontinuities, the numerical 
spread is likely to be larger than the physical one. 

Underestimates as well as overestimates of the radia- 
tive cooling at the contact discontinuities may occur when 
in the numerical code a cooling curve calculated in a regime 
of ionization equilibrium is implemented (as we actually do). 
A hot plasma undergoing strong radiativ e losses may cool 
at lo wer rat es j^u^herland fc Dopitall993F) or at higher rates 
(e.g. lBorkowskiet^lTn^9oir relative to the condition of ion- 
ization equilibrium, depending on its thermal history. 

In order to test the sensitivity of our models to the 
amount of radiative losses, we rerun model Draco-50 reduc- 
ing the radiative cooling by a factor of ten. Also in this case 
the great majority of the gas remains bounded to the galaxy. 
In particular ~ 70% of the gas remains inside the galactic 
volume, while only 30 % (with some oscillations) stays in the 
stellar region. Given the reduced effectiveness of the radia- 
tive cooling, the hot gas can expand more efficiently pushing 
most of the ISM at larger distances and creating a bubble 
encompassing a large fraction of the stellar volume. We thus 
conclude that the ability to retain the gas is a robust prop- 
erty of our model and depends essentially on the prolonged 
and thus low SNe II rate. As an argum ent in favour of this 
conclusion, we quote iMori et alJ {2002) who worked out a 
3D simulation of a subgalactic object with a potential well 
and an amount of ISM quite similar to our, with a gas bind- 
ing energy four times larger but a number of SNe II only two 
time higher. Despite this more disadvantageous proportion, 
most of the gas is blown away because the SNe II are made 
exploding all together instantaneously. 

7.3 Chemical evolution 

Although the majority of the metals expelled by SNe is still 
inside the dark matter halo or in its neighbourhood (and is 
still bounded to the galaxy), only a small fraction (~ 18%) 
is located inside the stellar region. This is the amount of 
metals which contributes to the metallicity of the form- 
ing stars. Thus, being the chemical elements observed in 



stars the main observational constraints in chemical evolu- 
tion models, assuming that the majority of the SN products 
stay in the observable zone may lead to infer non realis- 
tic galactic properties and/or stellar properties. No matter 
if metals are lost by g alactic winds (as assumed by some 
dSph s chemical models iRobertson et al.ll2005l : ICarigi et alJ 
2002) whether are pushed at higher distances, our finding 
is constitent with the assumption that the majority of the 
metals do not enrich the newly forming stars. 

At the early stage of the evolution the [Fe/H] distribu- 
tion in the ISM is rather spread. As the number of the bursts 
increases with time, the metals expelled by SNe become dis- 
tributed more and more homogeneously. This behaviour is 
reflected in the MDF which narrows around the mean value 
([Fe/H]}, which increases with time. At the end of the sim- 
ulations we obtain a mean value of ([Fe/H[) ~-1.65 with 
a spread of ~ 1.5 dex j Lehnert et alj|l99jj; ISh etrone et alJ 
l200ltlAparicio et aliEoOlt iBellazzini et al.M2002l) . while the 
maximum value of the distribution occurs at [Fe/H] =-1.5 
for the reference model Draco-50, and at [Fe/H]=-1.6 for 
Draco-25. 

We can also reproduce satisfactory the general ob- 
served behaviour of [O/Fe] in the [0/Fe]-[Fe/H] diagram. 
In dSphs this diagram shows a break in [O/Fe] occurring at 
[Fe/H]~ —1.6, with an approximate constant level of [O/Fe] 
at lower metallicities and a monotonic decrease at larger 
metallicities. In agreement with previous chemical evolution 
models (e.g. LM and IA) we regard the stars populating the 
plateau as polluted mainly by metals produced by SNe II. 
The shallow decline of the plateau at larger values of [Fe/H] 
is due to the contribution to the iron enrichment by SNe la 
which increases with time. 

Instead, we attribute the decreasing branch in the 
[0/Fe]-[Fe/H] diagram to the low value of the porosity of 
SNIa remnants: given the low SNIa rate, these remnants are 
located quite apart one from another, and the iron ejected 
by SNe la is distributed rather inhomogeneously through the 
stellar volume (cf. Fig. lMH . As a consequence, stars forming 
in the (relatively small) volume occupied by SNIa remnants 
have a ratio [O/Fe] lower than those forming elsewhere. This 
effect can naturally account for the stars with low [O /Fe] and 
high [Fe/H], without neces sarily invoking a co mplex SFH as 
suggested for some dSphs dCarigi et alJl2002T) . 

In our models the plateau and the descending branch 
become populated at the same time, although with different 
proportions. As the galaxy evolves, the break moves toward 
larger values of [Fe/H]. This interpretation of the [O/Fc]- 
[Fe/H] diagram differs from previous chemical models in 
which the break appears only after a specific time, without 
changing its position later on. LM interpret the break as a 
sign of the occurrence of a galactic wind. IA, who adopt a 
closed-box model, associate the break to a sudden switchover 
of the SNe la explosions at a specified time ti a ~ 1.5 Gyr. We 
claim that there is no direct link between the break presence 
(and position) and ti a . This conclusion is reinforced by the 
time evolution of the SNIa rate (cf. the solid line in Fig. 
which, after a rapid rise before ~ 0.5 Gyr, varies only by a 
factor of 2 over 2.5 Gyr; this smooth temporal profile can 
not justify the break in the [a/Fe]-[Fe/H] diagram (unless a 
drastic reduction in the star formation rate occurs). 

In the end, one possible discrimination between the LM 
scenario and the our would be given by the possible observa- 
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tion of some stars located "below" the plateau in the [a/Fe]- 
[Fe/H] diagram, i.e. showing a [O/Fe] ratio lower than that 
of the plateau at low values of [Fe/H]. 

We also stress t hat observations of "transition" galaxies 
jGrebel et al 1 12003) which have not yet lost all their gas 
content may discriminate between the LM scenario and the 
our. In the former one, infact, these galaxies are expected 
to show only the plateau in the [O/Fe] distribution, while 
in the latter one the descending branch would be already 
present at low values of [Fe/H]. 



7.4 Summary 

In conclusion, our major findings are the following: 

1) The long duration of the star formation attributed 
to Draco (and the other local dSphs too), together with the 
high effectiveness of the radiative losses, exclude that the 
galaxy got rid of its gas by an internal mechanism such as a 
galactic wind. The gas removal to complete star formation 
and to evolve to a gas-poor system should result from exter- 
nal mechanisms such as ram pressure stripping and/or tidal 
interaction with the Galaxy. 

2) Although the SN ejecta remain gravitationally 
bounded during the star formation, yet only a low fraction 
(~ 18%) stays in the region where star forms. This effect 
should be taken into account in chemical evolution models. 

3) Our models succeed in reproducing the [Fe/H] dis- 
tribution function of the stars. In agreement with observa- 
tions, we find a mean value ([Fe/H]) = —1.65 with a spread 
of ~ 1.5 dex. 

4) We can also satisfactory reproduce the observed 
[O/Fe] vs [Fe/H] diagram. 

5) Contrary to the usual interpretation, we rule out any 
particular relation between the [a/Fe] break position in the 
[«/Fe]-[Fe/H] plane and the onset of SN la explosions. Meth- 
ods intended to determine ti a in this way might be incorrect. 

6) In agreement with observations, our models de- 
velop initially moderate metallicity gradients which become 
weaker and weaker as the galaxy evolves. 

7) Finally, the above results are slightly dependent on 
the particular star formation history adopted. 
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APPENDIX A: 

A SNR expanding in a homogeneous uniform medium 
evolves initially following the Sedov solution during which 
radiative losses can be neglected. With time, however, such 
losses eventually cause the formation of a dense shell, and 
the high pressure of the hot interior pushes this shell (the 
snowplow phase). Assuming an explosion energy of 10 51 
erg, the time at wh ich the Sedov phase ends is given by 
ICioffi fc Shui] Jl99ll) t = 1.49 x 10 4 n~ 4/7 C~ 5/14 yr, where 
no is the density of the ambient medium and ( is the 
ISM metal abundance in solar units. The SNR stalls when 
its expansion velocity drops to the the local sound speed 
Co; the stalling radius is given by ICioffi fc Shulll (Il99lf) 
R Bt — 59rio 6/ ' 35 p 4 _1 '' 5 ^- 2 / 35 p C) i n terms of the pres- 
sure of the ambient medium P4 = Pj (&b X 10 4 cm -3 K) 
(where fee is the Boltzman constant), after a time t st = 
1.38ng- 33 P 1 T 7/10 C" 2/35 Myr. Eventually, the SNR cavity is 
refilled by the ISM on a timescale of the order of tsxi = 

-Rst/Co- 

It is known that a single SNR does not contribute sig- 
nificantly to the ISM energization; at the end of its evolu- 
tion only a few percent of the explosion energy has been fed 
into the s urrounding gas, while the rest has been radiated 
away (e.g. iBradamante et al"]ll998T) . However, when many 
SNe explode with a rate per unit volume S high enough, 
their combined effect may lead to the thermalization of a 
very high fraction (close to 100%) of the explosion energy. 
In fact, if many SNRs collide with each other the shells 
fragment and the hot gas of their interior merges forming 
a dilute medium. Repeated supernova explosions are espe- 
cially effective in reheating this medium because the SNR 
cooling time is relatively long in a low density region. A 
considerable fraction of gas may thus be driven out of the 
galaxy. To understand when this is effectively the case, one 
has to calculate the collision time tint- This time can be 
defined as the time required for SNRs forming at the rate 
5* to fill up the volume of space. If the expansion law of a 
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remnant is given by R s = At a (with A and a constants) , 
then fcf. lLarsonlll974l : lMelioli fc de Gouveia Dal Pinoll2004l) 
tint = (^-) 1/(3q+1) . Unless for very low values of no, the 
Sedov phase represents only a short stage in the evolution of 
a SNR whose bulk expansion can be described by the snow- 
plow law. In this case we have R B = 1.35n^ 13 ^ 49 C~ 2 ^ 49 tyr 7 
pc, where t yl is the time in years JCioffi fc ShuUlll99ll) . 

In order to make some quantitative consideration con- 
cerning our standard model, we assume no = 0.1 cm -3 , 
Pi — 0.1, 7?.snii = 3.7 x 10~ 6 yr _1 and a mean value 
S = 3.28 x 10" 15 pc -3 yr _1 . Then we have i int = 20.1C° 07 
Myr, which is longer than t st . This means that the SNRs do 
not interact one with another during their evolution because 
they stop their expansion before colliding. Despite this re- 
sult, with a sufficiently high rate S, and a long lifetime tmi, 
the SNR cavities would still necessary overlap, connect, and 
form a network of tunnels. One can characterize this topol- 
ogy in terms of the porosity parameter Q = SVsttmi, where 
V Bt = (47r/3)i? 3 t . The condition Q > 1 implies that SNRs 
overlap. With our assumptions we get Q = 0.5£ -8//35 . We 
thus conclude that in our model a combined action of SNRs 
pressurizing the ISM can be obtained for £ < 0.05. Assum- 
ing £ = 0.01 we calculate that SNRe II fill most of the stellar 
volume after t ~ 12 Myr. 



